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SURFACE  AND  INTERNAL  SIGNATURES  OF  ORGANIZED 
VORTEX  MOTIONS  IN  STRATIFIED  FLUID 

1.  Introduction 

Organized  vortices  form  in  nature.  The  formation  of  wake  vortices  in  the  stratified 
atmosphere  and  ocean  has  been  observed  (Chopra  &  Hubert  1965,  Lyons  &  Fujita  1968, 
Pao  &  Kao  1976,  Thomson  tt  al  1977),  analyzed  (e.g.,  Chopra  1973),  and  reproduced  in 
laboratories  (Barnett  1972,  Pao  Kao  1977).  The  atmosphere  vortices,  normally  identified 
by  the  shallow  stratocumulus  clouds  below  a  capping  inversion,  form  when  obstacles  (such 
as  island  mountains)  are  higher  than  the  inversion  with  a  large-scale  parallel  wind  in  the 
lower  atmosphere. 

Coherent  structures  in  stratified  fluids,  especially  solitary  waves  produced  by  the  col¬ 
lapse  of  a  mixed  region  of  fluid  (Wu  1969,  Kao  &  Pao  1979),  and  horizontal  vortices 
evolving  in  the  late  wakes  of  a  submerged  moving  object  (Pao  &  Kao  1977,  Lin  &:  Pao 
1979),  have  puzzled  many  researchers  because  of  the  striking  phenomena  these  flows  reveal 
and  their  distinct  characteristic  difference  from  those  in  homogenous  fluids.  The  puzzling 
stems  from  at  least  two  reasons.  First,  the  coherent  structural,  tnat  are  intrinsic  in  inter¬ 
nal  waves  or  wakes  exist  only  in  a  stratified  environment  but  not  in  a  homogeneous  one. 
Second,  flew  structures  in  internal  waves  or  wakes  dissipate  very  slowly.  Once  generated 
they  either  propagate  very  long  distance  or  persist  very  long  time. 

Density  inhomogeneity  and  gravity  play  a  very  subtle  role  in  the  dynamics  of  stratified 
fluid  motions.  Density  inhomogeneity  by  itself  has  very  minor  effects  on  flows.  Neither 
does  gravity  have  influence  on  fluid  flows  in  a  homogeneous  medium.  The  presence  of 
the  inhomogeneity,  however,  allows  gravity  to  play  a  role  in  the  dynamics  of  fluid  mo¬ 
tions.  Theoretically  speaking,  the  interaction  of  the  density  inhomogeneity  with  gravity 
or  with  other  force  fields  is  present  as  a  second  order  singularity  in  the  equations  for  flow 
stability,  and  therefore,  can  not  be  ignored  no  matter  how  small  the  density  inhomogene¬ 
ity  is  (Mile  1961,  Howard  1961,  Fung  &  Kurzweg  1975).  Such  interaction  produces  an 
oscillation  mechanism  that  dominates  the  motion  of  the  coherent  structures  in  an  inhomo¬ 
geneous  environment.  This  the  why  the  very  weakly  stratified  atmosphere  and  ocean  can 
support  internal  waves,  wakes,  and  otiiei  organ  Led  duv.  atiuciurcs  that  do  not  exist  in  a 
homogeneous  environment. 

Muuucript  approved  April  28,  1993. 
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Vortex  motions  in  a  gravity-stratified  fluid  are  further  complicated  by  three  other 
factors;  the  centrifugal  force  field  generated  by  the  rotation  of  the  flow,  the  radial  density 
gradient  induced  by  the  pressure  constraint,  and  the  interaction  of  the  two  force  fields  with 
the  density  inhomogeneities  in  both  the  gravitational  and  centrifugal  directions. 

The  rotation  of  vortices  plays  a  dual  role  in  the  stability  of  the  organized  flow  struc¬ 
tures.  While  the  angular  velocity  gradient  produces  shear  effects  which  always  destabilize 
the  flow,  the  angular  velocity  itself,  interacting  with  the  density  inhomogeneity,  creates  a 
centrifugal  force  field  that  can  either  stabilize  or  destabilize  the  flow  field  (Fung  1983). 

As  required  by  the  balance  of  pressure  in  all  directions,  the  angular  velocity  of  vortex 
motions  automatically  induces  density  gradients  in  the  radial  direction.  This  induced 
density  gradient  allows  the  centrifugal  force  to  play  a  role  in  the  dynamics  of  the  flow 
structures  similar  to  what  the  gravity-stratified  density  gradient  does  in  the  gravitational 
force  field,  and  becomes  another  driving  force  for  the  motion  of  the  coherent  structures  in 
a  stratified  environment. 

As  a  result  of  the  above  interactions  of  the  force  fields  and  density  inhomogeneities  in 
more  than  one  direction,  any  organized  flow  structures,  if  existing,  will  have  to  satisfy  all 
the  required  conditions.  This  includes  the  stability  condition  in  the  gravitational  direction, 
the  stability  condition  in  the  centrifugal  direction,  and  the  pressure  balance  condition 
imposed  by  the  pressure  constraint  (Fung  1986,  1991). 

The  presence  of  free  surface  in  vortex  motions  in  a  stratified  fluid  even  further  com¬ 
plicates  the  already  complex  dynamic  system.  The  barotropic  mode  prompted  by  the  free 
surface  is  fast  moving  and  must  be  treated  numerically  for  computational  stability.  To 
efficiently  handle  this  fast  moving  inode,  a  split-explicit  numerical  scheme  (Chang  19S4, 
1985)  was  used  to  separately  integrate  different  modes  in  the  ocean.  In  addition  to  this 
barotropic  mode,  the  change  of  the  surface  height  also  complicates  the  balance  of  the 
centrifugal  and  gravitational  forces  and  thus  their  interaction  with  density  distribution. 
Besides,  the  perturbation  of  the  free  surface  Flight  is  difficult  to  observe  and  measure  both 
in  laboratories  and  over  the  ocean,  and  thus  further  hampers  the  understanding  of  the 
phenomenon. 
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As  observed  by  Pao  &  Kao  (1977)  in  their  laboratory  experiment,  when  a  sphere  is 
towed  through  a  stratified  fluid,  vortices  are  first  shed  without  any  preference  in  direction. 
The  region  immediately  behind  the  towed  body  is  in  a  turbulent  state  with  a  considerable 
amount  of  mixing.  This  state  of  turbulence  is  quickly  inhibited  by  the  density  stratification 
in  the  vertical  direction.  The  mean  motion  in  the  horizontal  direction  then  dominates  the 
wake.  The  flow  first  meanders  and  then  gradually  rolls  up  into  vortices.  The  resultant 
pattern  formed  in  the  late  wake,  even  though  inherently  three-dimensional,  is  confined 
within  a  relatively  thin  horizontal  layer,  and  is  reminiscent  of  the  two  dimensional  Karman 
vortex  street  if  observed  from  the  gravitational  direction. 

Most  of  the  modelling  of  wake  vortices  are  conducted  with  a  two-dimensional  model 
(e.g.,  K&miadakis  Triantafyllou  1989)  in  which  the  effects  of  stratification  and  hori¬ 
zontal  mass  divergence  are  not  properly  considered.  The  wake  vortices  generated  in  such 
a  model,  though  appear  like  the  Kaman  vortex  streets,  may  have  a  completely  different 
mechanism  of  formation.  As  will  be  discussed  in  Section  5  of  this  paper,  no  corresponding 
two-dimensional  model  can  adequately  describe  the  horizontal  vortices  in  stratified  fluids 
because  the  stratification  can  not  be  properly  parameterized. 

This  paper  will  (1)  briefly  describes  a  three-dimensional  numerical  model  with  a  free 
surface  that  is  capable  of  generating  the  vortex  motion  behind  a  submerged  moving  body 
or  other  organized  flow  structures  in  a  stratified  environment;  and  (2)  propose,  based  on 
the  numerical  results  obtained,  a  mechanism  for  the  generation  and  growth  of  horizontal 
vortices  in  stratified  fluids.  It  is  not  our  intention  to  realistically  simulate  the  detailed 
flow  structure  near  a  submerged  moving  body.  Instead,  our  interest  is  the  development  of 
organized  flow  structures  after  the  initial  wake  turbulence  induced  by  submerged  moving 
objects  has  dissipated.  We  have  limited  the  scope  of  this  paper  to  the  generation  and 
growth  of  organized  flow  structures  in  the  stratified  late  wake.  Specifically.  Section  2 
will  present  the  formulation  of  model  dynamics,  turbulent  mixing,  mode  splitting,  and 
numerical  techniques.  Section  3  will  describe  the  near  field  velocity  profile  to  generate 
the  horizontal  vortices.  Section  4  will  discuss  the  p  minary  results  of  such  vortices  and 
examine  the  relationship  between  the  flow  structure  and  the  induced  free  surface  signature. 
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Section  5  will  present  a  mechanism  for  the  generation  and  evolution  of  horizontal  vortices 
based  on  the  numerical  results.  Section  6  will  conclude  our  findings. 


2.  The  3D  Numerical  Model 

A  three-dimensional  time- dependent  nonlinear  model  is  constructed  to  study  the  wake 
vortices  in  a  stratified  fluid.  The  flow  is  assumed  to  be  incompressible  and  hydrostatic.  The 
Boussinesq  approximation  is  applied.  The  governing  equations,  the  turbulent  closures,  the 
splitting  of  barotropic  and  baroclinic  modes,  the  grid  structure,  the  boundary  and  initial 
conditions  of  the  model  will  be  discussed  in  this  section.  The  model  describes  the  motion 
of  a  stratified  fluid  in  a  tank-like  confine  in  the  Cartesian  coordinates  with  a  flat  rigid 
bottom,  a  free  surface,  and  close  or  open  side  walls.  The  x-direction  is  the  direction  of 
towing,  the  y-direction  is  cross-stream,  and  the  2-direction  is  vertical  and  opposite  to  the 
gravity.  As  shown  by  Figure  1,  the  fluid  in  the  model  has  an  average  height  of  H  with  a 
free  surface  perturbation  h  such  that  h  <<  H. 

a.  Governing  Equations 

The  horizontal  equations  of  motion  in  the  r-  and  y-directions  of  the  model  are 


du  du  du  du 

¥  +  ua;  +  ''^  +  u'aI 


1  dp 

—  J:  +  Fh(u)  +  F'(v) 
po  ox 


(1) 


dv  dv  dv  dv 
dt  +  Udx+Vdy+Wdz 


-~  +  Fh(v)  +  Ft(v) 
po  Oy 


(2) 


Here  u,v,  and  w  are  velocity  components  in  the  i,y,  and  ;  directions,  respectively,  and 
po  is  a  density  constant.  Functions  F/,  and  F,  represent  the  effects  of  turbulent  mixing 
from  the  horizontal  and  vertical  fluxes.  The  equation  of  motion  in  r-uirection.  under  the 
hydrostatic  assumption,  is  reduced  to 


dp 

d : 


=  -PQ 


(3) 
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where  p  is  the  density  of  the  fluid.  For  the  far-field  simulation  as  in  the  present  case, 
organized  vortices  form  long  after  the  initial  disturbances  generated  by  the  moving  object 
have  dissipated.  The  hydrostatic  assumption  is  justified  at  this  stage  when  the  buoyancy 
effect  due  to  density  stratification  dominates  the  vertical  motion.  The  mass  conservation 
of  the  model  is 


+  +  =  Fh{p)  +  FM 


w 


For  incompressible  fluids,  the  continuity  equation  is 


du  dv 
dx  +  dy  + 


dw 


=  0 


(5) 


With  the  rigid  bottom  where  the  vertical  velocity  vanishes,  i.e.,  w  =  0  at  z  =  0,  an  vertical 
integration  of  (5)  from  the  bottom  to  the  surface  of  the  fluid  (2  =  H  4-  h)  yields 


l 


H+h  dw 

-r-dz  =  w(H  +  h)  -  u>(0)  =  w{H  4-  h) 
Oz 


-l 


dx  dy' 


dx 


,  d{v)  du,  dv, 

+  'aT)  +  ',(a7  + W* 


where  h  is  the  surface  height  perturbation  and  the  depth  average  (d)  of  any  function  0  is 
defined  as 


where 


d  =  d'  +  (o) 


(6) 


The  prognostic  equation  for  the  surface  height  perturbation  h,  after  rearranging,  is  ob¬ 
tained  as 


o 


dh  tT,dv  .  dv  r  d  d  i 

af  =  +  '  la7<hu*)  +  aJ(H 


(8) 


where  u,  and  v ,  are  surface  current  velocities  with  negligible  variations  across  A. 

&.  Turbulence  Mixing 

There  are  many  ways  to  model  the  turbulent  mixing.  The  available  parameteriza- 
tions  vary  horn  the  simple  /C-theory  with  a  constant  diffusion  coefficient  to  higher  order 
turbulent  closures.  Direct  numerical  simulations  require  no  parameterization  but  demand 
much  more  time  and  capability  in  computing.  To  simplify  the  computation,  we  opt  for  a 
somewhat  conservative  turbulent  closure  just  enough  for  the  model  to  form  vortices. 

For  the  horizontal  diffusion,  Smagorinsky’s  formula  is  used  to  accept  contribution  from 
the  rate  of  deformation  of  fluid  elements  as  follows  : 


Here  k  is  the  von  Karman  constant,  Ax  and  Ay  are  grid  sizes,  and  the  deformation  fields 
are  defined  as 


G 


Kv  -  (A z)2 \/[(^)J+(^)5](l-i?,),  for  1  >  R,  >  0. 


(12) 


where  the  Richardson  number  is  defined  as 


Rt  = 


and  the  Brunt-Vaisala  frequency  is  defined  as 


N2  = 


9  dp 

p  d: ' 


(13) 


The  vertical  mixing  vanishes  as  the  gradient  Richardson  number  equals  to  or  larger  than 
unity.  The  vertical  turbulent  mixing  function  in  (1)  and  (2)  then  defined  as 


Ft{u,v,p) 


d:  “  d:  J  ’ 


(14) 


c.  Splitting  of  lnttmal  and  External  Modes 

The  governing  equations  (1)  to  (5)  and  (8)  contain  the  internal  gravity  (baroclinic) 
and  external  gravity  (barotropic)  modes.  For  a  certain  depth  and  density  stratification  in 
the  ocean  or  in  wave  tanks,  the  characteristic  phase  speed  of  the  externa]  gravity  mode 
c  =  vW  is  about  one  or  two  orders  of  magnitude  larger  than  the  first  internal  gravity 
mode.  In  various  configurations  of  our  model,  the  average  depth  H  ranges  from  30  to 
70  cm,  giving  a  phase  speed  of  171  to  262  cm/ sec  for  the  external  modes.  The  maximum 
internal  phase  speed,  however,  is  dictated  by  the  towing  speed,  which  is  4  cm/ sec.  It 
is  expedient  to  separate  the  governing  equations  into  a  barotropic  part  and  a  baroclinic 
part  and  to  numerically  integrate  them  at  different  time  intervals  to  achieve  computational 
efficiency  (Chang  1984,  1985). 

Now  we  rewrite  (1)  and  (2)  as 


du_  1  dp 

dt  ‘  p0  di 


(15) 


( 


(16) 


where 


—  =  5  --  JL^£ 
Of  po  Oy 


Ou  Ou 

uoI  ~v!h 


+  F*(u)  +  F,(u) 


Ov  Oi>  Ov  ,,  _  .  . 

B  =  ■“&  -  '%  -  u,57  +  n<”> + 

For  any  point  in  the  fluid  at  the  depth  of  H  +  h  —  2,  the  total  pressure  is  equal  to 

/H  +  /»  ?H  +  h 

pdz'  =  g  J  pdz'  +  gj"  pdz'  (17) 

The  value  of  p(z)  in  (17)  contains  a  large  portion  of  static  pressure  p0g(H  —  z)  which  does 
not  contribute  to  the  dynamics  of  the  motion.  Defining  a  reference  pressure 


pr{z)  =  p{z)  -  p„g{H  -  z) 
and  substituting  into  (17),  we  have 


P  r 


H 


e'dz'  +  g(t)(H  -  r)  +  ( p ,  ~  p0)gh , 


(IS) 


where  pt  is  the  fluid  density  at  the  free  surface,  e  =  p  —  p0  the  density  deviation,  (c)  the 
depth  average  of  e  as  defined  in  (6),  and  t'  —  t  —  (e)  as  in  (7).  To  arrive  in  (18)  we  have 
assumed  that  the  density  variation  across  h  is  negligible. 

The  second  term  in  (18)  is  the  static  pressure  produced  by  the  depth  average  density 
variations  in  stratified  fluids.  The  last  term  in  (18)  contains  the  contribution  from  the 
perturbation  of  the  free  surface,  a  depth-independent  pressure  force  exerting  on  th^  fluid. 
Depth  averaging  (18)  yields  the  depth  average  reference  pressure  and  its  deviation  as  : 


S 


(19) 


(pr)  =  r {<)9h  +  s(«')  +  ( p .  -  p0)gh , 


and 

Pr  =  Pr  -  (Pr). 


(20) 


The  depth  average  pressure  gradient  force,  V(pr),  exerts  uniformly  throughout  the 
fluid  column  and  excites  the  bnrotropic  mode  in  the  fluid.  The  deviation  of  the  pressure 
gradient  force,  on  the  other  hand,  excites  only  the  baroclinic  mode  and  contributes  more  to 
the  detailed  internal  structures  of  the  organized  vortex  motion.  Taking  the  depth  average 
of  (11)  and  (12),  we  separate  the  prognostic  equations  (1)  and  (2)  into  the  barotropic  set 


d(“)  _  ,  A\  _  1  d(pr) 

dt  ~  {  }  Po  dx 


(21) 


W  __  /  D\  1  dWr) 

sr  - {B)  -  7.~sr- 

and  the  baroclinic  set 


5t/  _  ,  1  dp'r 

dt  “  Po  dx  ’ 


(23) 


dt  ~  °  ~  Po  dy  ’ 


(24) 


dt 


d('  dt' 
Udx  V'dy 


+  f A  (  «  )  +  Fi :(0- 


(25) 


One  can  of  course  further  separate  the  internal  modes  such  as  in  Madala  k  Piacsek  (1977) 
or  in  Gill  (19S4).  however,  the  small  differences  in  phase  speeds  among  baroclinic  inodes 
yields  a  diminishing  computational  benefit  for  further  splitting. 
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d.  Grid  Structure  and  Finite  Differences 

Figure  2  illustrates  the  three-dimensional  grid  with  variables  staggered  in  Arakawa  C 
configuration  (Arakawa  &  Lamb  1977).  In  the  C  configuration,  mass  variables  such  as  p 
and  p  are  defined  at  the  full  points  (•  points).  The  u-components  (o  points),  v-components 
(o  points),  and  u>-components  (A  points)  are  defined  at  a  half  of  a  grid  length  offset  from 
the  mass  points  in  the  x-,  y-,  and  ^-directions,  respectively.  For  the  results  to  be  reported 
here,  a  grid  with  163  x  41  x  15  mass  points  in  the  x-,  y-,  and  r-directions  is  used  with 
a  spatial  resolution  Ax  =  Ay  =  A z  =  2  cm  in  the  unstretched  region.  The  grid  can  be 
stretched  with  various  ratios  in  all  three  directions.  In  the  present  simulation,  the  grid 
is  stretched  only  in  the  x-direction  for  x  >  280  cm  with  a  ratio  equal  to  1.3  :  1.  This 
stretching  allows  us  to  simulate  the  slowly  evolving  vortices  while  keeping  the  submerged 
object  moving  downstream,  and  to  avoid  the  possible  adverse  effect  from  the  downstream 
reflecting  gravity  waves  on  the  wake  vortices.  The  spatial  differences  are  of  second  order 
accuracy. 

e.  Boundary  and  Initial  Conditions 

A  free- slip  boundary  condition  is  used  for  the  horizontal  velocities  at  the  rigid  bottom 
where  the  vertical  velocity  w  =  0.  The  free-slip  condition  is  justified  since  our  region 
of  interest  is  away  from  the  bottom  boundary.  At  the  two  walls  along  the  direction  of 
the  moving  object,  a  Neumann  boundary  condition  is  used.  To  avoid  the  reflection  of 
hump  waves  that  travel  with  the  vehicle  from  the  downstream  wall,  the  stretched  grid 
accompanying  with  larger  Ks  acts  as  an  effective  absorbing  region.  The  initial  density 
distribution  is  statically  stratified  with  a  uniform  strength  N  =  0.01/sec.  We  have  also 
tested  a  case  for  homogeneous  fluids  ( N  =  0)  and  produced  no  organized  flow  patterns. 

/  Temporal  Integration 

Following  Chang  (1984,  1985)  in  his  split-explicit  leapfrog  scheme  for  temporal  in¬ 
tegration,  wr  max :h  the  barotropic  and  baroclinic  sets  of  equations  forward  with  a  time 
interval  of  0.003125  and  0.25  sec,  respectively.  The  computation  cycle  starts  with  calcu¬ 
lating  the  hydrostatic  reference  pressure  and  vertical  velocity  u-  based  on  the  current  state 
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of  the  flow.  We  then  evaluate  the  transport  and  diffusion  terms  A  and  B  in  (15)  and  (16). 
and  obtain  (A)  and  ( B )  by  vertical  integration.  The  barotropic  set  set  of  equations  is 
subsequently  integrated  80  times  with  the  depth  average  pressure  gradient  force  updated 
each  time  step  by  a  new  h  in  (19),  (21),  and  (22).  The  values  of  A'  and  B'  remain  constant 
during  these  80  steps.  The  baroclinic  set  of  equations  are  then  integrated  with  the  big 
time  interval.  At  last  we  complete  one  computational  cycle  by  updating  all  the  dependent 
variables. 

g.  Tracer  Injection 

We  also  introduce  a  transport  equation  for  an  inert  tracer  with  concentration  Q  and 
a  source  term  q  as  follows: 


The  purpose  of  this  equation  is  to  reveal  the  evolution  of  possible  organized  flow  structures 
in  the  stratified  fluid. 


3.  Description  of  the  near  field  velocity 

Modelling  the  near-field  flow  generated  by  a  moving  object  is  extremely  difficult.  The 
complex  boundary  conditions  inherent  in  the  object  geometry  and  the  high  resolution 
required  for  the  thin  boundary  layer  of  the  object  demand  great  computing  resources. 
Substituting  the  object  with  a  corresponding  flow  field  within  the  computational  domain 
therefore  becomes  a  very  attractive  approach  to  circumvent  these  difficulties.  Abernathy 
&  Kronauer  (1962)  first  applied  this  idea  to  simulate  the  instability  of  two  parallel  vortex 
sheets.  Recently  Triantafyllou  &  Karniadakis  (1990)  used  this  approach  to  simulate  the 
vortex  streets  by  the  spectral  element  method.  We  follow  the  same  strategy  by  producing 
a  flow  field  that  is  at  a  finite  distance  from  the  immediate  boundary  of  a  towed  sphere. 
This  approach  fits  in  our  region  of  interest  in  which  the  organized  flow  structures  develop 
in  a  region  far  away  from  the  moving  object. 
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To  emulate  the  towed  body  as  in  several  laboratory  experiments  (Lin  £:  Pao  1979, 
Pao  &  Kao  1977),  a  spherical  region  with  a  diameter  D  =  6  cm  is  prescribed  in  the  model 
to  represent  a  towed  object.  Within  the  region,  the  flow  velocity  is  artificially  relaxed 
towards  a  target  velocity.  As  shown  in  Figure  3,  the  target  velocity  in  the  x-  direction 
(u0)  is  equal  to  that  of  the  prescribed  moving  speed  of  the  spherical  region  (or  the  object), 
4  cm/sec.  The  target  velocity  in  the  y-  direction  (v0),  flowing  radially  inwards  to  simulate 
a  convergent  flow  field,  is  equal  to  a  quarter  of  the  moving  speed  within  the  spherical 
region.  The  flow  at  those  grid  points  within  the  radius  of  the  sphere  is  accelerated  toward 
the  target  velocity  by  adding  a  forcing  term  to  the  horizontal  equations  of  maion  (15)  and 
(16)  as  follows  : 


du  _  1  dp 

di  ~  A~  JoTt 


-  A(u  -  u0) 


(27) 


1  dp 

po  dy 


-  A(i>  -  v0) 


(28) 


Here  A  is  the  e-folding  time  scale  for  the  forcing  and  is  set  to  1/sec  in  the  present  example. 
The  vertical  velocity  w  is  computed  by  the  continuity  equation  (5)  which  produces  a 
vertical  convergent  flow  field  to  simulate  the  vertical  wake  collapse  phenomena. 

The  movement  of  the  spherical  region,  with  a  speed  of  4 cm /sec  in  the  downstream 
x-direction,  is  described  by  the  trajectory  of  its  center  located  right  in  the  middle  of 
the  y-z  plane.  To  simulate  the  unsteadiness  of  the  moving  body,  we  further  introduce 
a  random  perturbation  in  the  y-component  of  the  center  of  the  sphere  with  a  maximum 
perturbation  amplitude  less  than  0.2  %  of  the  diameter  of  the  spherical  region.  The 
introduction  of  the  random  perturbation  is  realistic  because  of  the  random  background 
noise  and  vibrations  in  our  environment.  An  amplitude  of  0.001  cm  is  much  smaller 
that  what  can  be  achieved  in  laboratories.  This  near  field  acceleration  continues  to  be 
enforced  as  the  prescribed  sphere  travels  downs  ream  to  simulate  the  continuous  moving 
of  a  submerged  object.  7  urthermorc.  the  inert  tracer  is  released  ronstantly  from  inside  of 
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the  sphere  at  a  prescribed  rate  to  illustrate  the  evolution  of  organized  vortices  in  the  flow 
field. 


4.  Vortex  Evolution 

At  t  =  0  the  prescribed  sphere  starts  moving  at  x  =  7  cm.  To  illustrate  the  results, 
we  select  a  horizontal  plotting  window  that  covers  a  region  from  x  =  140  to  240  cm  and 
y  =  16  to  64  cm.  The  sphere  enters  from  the  left  of  the  window  at  t  ~  29  sec  and  exits  to 
the  right  at  t  ~  54  sec. 

The  x-y  plane  in  which  the  results  are  displayed  in  Figures  4  to  7  is  in  the  middle  of 
the  channel  ( z  =  15  cm)  and  contains  the  trajectory  of  the  center  of  the  moving  sphere. 
Figure  4  shows  the  evolution  of  the  current  vectors  in  the  middle  layer  at  different  time. 
At  t  <  120  sec,  the  horizontal  velocity  is  more  or  less  parallel  to  the  x-direction  as  a 
result  of  the  downstream  movement  of  the  spherical  body.  Later  at  t  =  360  and  600  sec, 
the  flow  starts  to  meander,  then  gradually  rolls  up,  and  finally  develops  it  to  organized 
vortices  with  alternating  patterns  as  shown  at  later  time  in  Figure  4.  The  corresponding 
maximum  vectors  for  the  horizontal  velocity  within  the  plotting  window  slowly  decreases 
from  0.291  cm/sec  at  f  =  120  sec  to  0.0835  cm  I  sec  at  i  =  1800  sec,  indicating  the  slow 
rate  at  which  the  organized  flow  structures  dissipate  in  the  stratified  fluid. 

The  tracer  concentrations  in  Figure  5  illustrate  the  evolution  of  the  coherent  vortex 
patterns  at  the  corresponding  time  as  in  Figure  4.  In  a  very  slow  evolution,  the  patterns 
meander,  roll  up,  and  finally  develop  into  vortices  in  the  same  manner  as  the  velocity 
does.  The  evolution  of  the  tracer  pattern  resembles  most  to  what  has  been  observed  in 
laboratory  experiments. 

The  density  anomaly  depicts  similar  behaviors  in  the  development  of  the  wake  vortices 
as  those  in  the  velocity  and  tracer  fields.  The  dash  lines  in  Figure  6  represent  negative 
density  deviations  (lighter  fluid)  from  the  initial  distribution.  It  is  apparent  that  the 
regions  of  lighter  density  are  located  near  centers  of  vortices.  This  is  consistent  with  the 
pressure  balance  requirement  for  stratified  vortex  motions  in  which  lighter  fluid  particles 
with  faster  velocity  can  embed  in  heavier  fluids  (Fung  1991). 
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Another  response  of  the  fluid  to  balance  the  centrifugal  force  is  the  surface  height 
perturbation.  The  surface  signatures  in  Figure  7  reflect  the  evolution  of  the  organized 
flow  structure  underneath  once  the  large  amplitude  hump  wave  disperses  away.  The  first 
panel  in  Figure  7  is  the  signature  at  t  =  120  sec  with  a  displacement  ranging  from  —4.5 
to  5  x  10~5cm.  The  surface  height  perturbation  at  this  time  reflects  only  the  gravity 
waves  due  to  the  passage  of  the  sphere.  As  also  illustrated  at  the  corresponding  time 
for  the  tracer,  velocity,  and  density  fields  in  Figures  4,  5,  and  6,  the  surface  signature 
reveals  no  organized  structures  within  the  flow  field.  As  the  initial  gravity  waves  dissipate 
and  the  vortex  patterns  start  evolving  in  those  fields,  however,  the  internal  flow  behavior 
leaves  very  discernible  surface  signatures  as  evident  in  the  free  surface  perturbation.  The 
corresponding  displacements  for  the  surface  motion  are  ranging,  respectively,  from  —2.8  to 
0.6  x  10-6 cm  at  t  =  360  sec,  to  —1.1  to  0.25  x  10~6cm  at  t  =  1S00  sec.  A  perspective  view 
for  the  free  surface  is  given  in  Figure  8,  showing  the  surface  perturbation  corresponding 
to  the  horizontal  vortices.  Except  the  one  at  t  =  120  sec  being  the  remnants  from  the 
passage  of  the  sphere,  the  induced  surface  signatures  respond  closely  to  the  evolution  of 
the  horizontal  vortices  underneath. 

Even  though  the  surface  signatures  shown  in  Figures  7  and  8  reveal  the  existence  of 
the  organized  flow  structures  under  the  free  surface,  they  may  be  affected  by  the  depth  of 
water  and  the  depth  at  which  the  submerged  object  moves.  In  other  words,  similar  internal 
vortex  structures  may  produce  different  surface  signatures  when  the  depth  of  w-ater  or  the 
depth  of  the  moving  object  vary.  These  behavior  will  be  discussed  in  a  separate  paper. 

Figure  9  shows  the  velocity  profile  of  a  vertical  plotting  window  that  covers  a  region 
from  i  =  140  to  240  cm  with  the  full  water  depth.  The  window  is  in  the  x-z  plane 
located  in  the  middle  of  the  computation  region  (y  =  40  cm)  which  contains  the  trajectory 
of  the  center  of  the  moving  object.  Contrast  to  the  organized  patterns  developing  in 
the  horizontal  plane,  little  evidence  of  organized  vertical  structures  is  shown  as  the  time 
evolves.  This  behavior  is  consistent  with  what  have  been  observed  in  experiments  for  the 
horizontal  vortices. 

The  density  distribution  at  t  =  120  sec  in  Figure  10  shows  the  density  profile  before 
the  horizontal  vortices  develop.  Layer  density  distributions  can  clearly  be  seen  at  this  stage 
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of  early  wake  where  turbulent  mixing  due  to  the  passage  of  the  sphere  is  quickly  suppressed 
by  gravity.  The  horizontal  vortices  then  start  to  evolve,  and  the  density  distribution  at 
later  time  responds  accordingly  to  balance  of  the  centrifugal  and  gravitational  force  fields. 
The  results  for  the  density  in  the  1-2  plane  show  a  different  facet  of  the  flow  structure 
corresponding  to  the  development  of  the  horizontal  vortices.  In  a  deeper  channel,  the 
vertical  flow  structure  will  further  reveal  the  coherence  of  the  organized  vortex  motion. 
Those  results  will  be  presented  in  an  upcoming  paper. 

5.  A  Proposal  Mechanism 

The  resemblance  of  the  Karman  vortex  streets  behind  a  bluff  body  and  the  horizontal 
vortices  in  stratified  wakes  would  lead  researchers  to  relate  one  to  the  other.  To  under¬ 
stand  the  generation  mechanism  of  the  horizontal  vortices,  one  would  intuitively  look  for 
a  mechanism  similar  to  that  of  the  shedding  for  the  Karman  vortex  streets.  The  mech¬ 
anism  extended  from  such  an  intuition,  however,  will  have  to  face  the  challenge  of  two 
fundamental  questions. 

First,  while  the  velocity  deficit  behind  a  bluff  body  is  responsible  for  the  shedding  of 
the  Karman  vortex  streets,  the  vortices  immediately  roll  up  behind  the  bluff  body,  become 
fully  developed  within  a  couple  diameters  behind  the  body,  and  rapidly  decay  as  they  are 
convected  downstream.  Large-scale  flow  structures  may  appear  as  the  secondary  vortex 
“street”  in  far  wakes  (Cimbala  et  al.  1988),  however,  they  are  not  steady  phenomena 
as  compared  to  those  vortices  in  the  primary  vortex  streets  or  in  the  horizontal  vortices. 
Even  if  appearing,  they  axe  only  characterized  by  lower  frequencies  as  the  continuing 
development  of  the  organized  structures  from  the  primary  streets.  If  a  similar  mechanism 
is  responsible  for  generating  the  horizontal  vortices  in  stratified  wakes,  it  can  not  explain 
why  the  horizontal  vortices  only  appear  long  after  the  initial  shedding  and  disturbances 
created  by  a  submerged  moving  object  have  dissipated.  That  is  a  time  long  after  the 
Karman  vortex  streets  would  decay  behind  a  bluff  body.  Also  different  from  those  of 
the  secondary  vortex  street,  the  flow  structures  of  the  horizontal  vortices  in  stratified 
fluids  always  evolve  steadily  from  small-scale  to  large-scale,  not  only  dominated  by  low 
frequencies. 
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Second,  Karman  vortex  streets  exist  in  both  homogeneous  and  stratified  fluids.  The 
density  inhomogeneity  may  have  influence  on  the  detailed  structures  of  the  streets  but  has 
no  bearing  on  their  generation  or  existence.  The  horizontal  vortices,  on  the  other  hand, 
exist  only  in  stratified  fluids.  The  shedding  mechanism  for  the  Karman  vortex  streets 
cannot  explain  why  the  horizontal  vortices  do  not  exist  in  homogeneous  fluids. 

In  view  of  the  two  questions  just  raised,  the  generation  mechanism  for  the  horizontal 
vortices  must  be  different  from  that  for  the  Karman  vortex  streets.  Instead  of  searching 
for  the  resemblance  in  shedding  mechanism  between  the  Karman  vortex  streets  and  the 
horizontal  vortices,  we  therefore  focus  on  the  fact  that  the  horizontal  vortices  develop  long 
after  the  initial  induced  turbulence  stage  is  over,  and  that  they  exist  only  in  stratified  fluids. 
From  our  numerical  results  obtained  in  the  present  simulation,  we  propose  a  mechanism 
in  which  the  density  plays  a  crucial  role  in  generating  the  horizontal  vortices  in  stratified 
wakes  as  follows. 

As  the  initial  disturbances  caused  by  the  passage  of  a  submerged  object  gradually 
dissipate,  the  buoyancy  effect  begins  dominating  the  vertical  motion  in  a  stratified  en¬ 
vironment.  The  hydrostatic  equation  referenced  to  a  fluid  element  with  density  pa  and 
pressure  p0  at  a  depth  of  z0  can  be  written  as 

~  Po)  =  (P  ~  Po)g  (29) 

where  p  —  pa  is  the  density  anomaly  which  interacts  with  the  gravity  to  produce  an 
buoyancy-induced  oscillation  with  respect  to  an  equilibrium  position  z0.  The  fluid  ele¬ 
ment  is  decelerated  by  the  gravity  above  the  equilibrium  position  and  accelerated  by  the 
buoyancy  below  the  equilibrium  position  as  described  by  : 


vu/  . 

—  <0  for  z  >  z0 

az 


—  >0  for  z  <  z0  (30) 

az 

The  situation  is  much  more  complicated  for  vortex  motions  in  which  a  centrifugal  force 
field  is  present.  Rotating  fluids  with  lighter  density  can  conceivably  be  embedded  in  an 
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environment  of  heavier  density,  however,  the  buoyancy-induced  restoring  force  remains  a 
dominating  factor  for  the  oscillation  in  stratified  late  wakes. 

The  transport  equation  for  the  7-component  vorticity  under  the  Boussinesq  assump¬ 
tion  reads 


dut  du>z  dw  dw  dw 

_+u_+„_+u,__„,_.+Wy_+Ui_+4, 


(31) 


where  and  w,  are  vorticity  in  the  x,y,  and  7-direction,  respectively,  and  $  repre¬ 

sents  turbulence  dissipation.  The  first  three  terms  on  the  right  hand  side  of  (31)  are  the 
stretching  of  the  three  vorticity  components,  a  vorticity  forcing  term  in  the  7-direction. 
The  last  term  in  (31)  is  the  result  of  the  turbulent  stresses,  a  diffusion  term  that  influences 
the  vorticity  distribution. 

If  the  vortex  lines  align  more  or  less  with  the  gravity,  the  vorticity  in  the  7-direction 
dominates  the  vorticity  distribution,  i.e.,  u,  >>  u>z,uy,  and  the  vorticity  equation  in  the 
7-direction  can  be  approximated  by 


du,'z  dw,  duj,  dw,  dw 

-ar  +  “ar  +  t'nr  +  u,arsu‘a7  +  * 


(32) 


The  vorticity  stretching  in  (32),  the  product  of  the  vertical  velocity  gradient  (VVG)  and 
the  vorticity  in  the  7-direction,  becomes  the  main  forcing  term  for  the  vorticity  tendency. 
For  a  positive  vorticity  distribution,  the  vorticity  increases  for  an  upward  accelerated  fluid 
element  and  decreases  for  a  downward  decelerated  fluid  element  according  to  (30).  It  is  vice 
versa  for  a  negative  vorticity  distribution.  We  will  use  our  numerical  results  to  illustrate 
the  effect  of  this  buoyancy-induced  oscillation  on  the  evolution  of  the  horizontal  vortices. 
The  vorticity  information  will  be  extracted  from  a  x-y  plane  in  the  middle  of  the  channel 
where  (32)  is  valid. 

When  the  submerged  body  passes  the  stratified  fluid  as  in  our  numerical  simulation 
and  in  laboratory  experiments,  a  three-dimensional  turbulent  wake  column  is  produced 
As  the  initial  turbulence  generated  by  the  moving  body  is  quickly  damped  by  stratification, 
this  wake  column  is  under  the  influence  of  the  gravity  or  buoyancy  force  in  the  vertical 
direction.  The  fluid  is  being  vertically  accelerated  or  decelerated  across  the  vake  column. 
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This  buoyancy-induced  acceleration  interacts  with  the  vorticity  within  the  wake  column 
in  a  way  described  by  equations  (30)  and  (32),  and  organizes  the  coherent  flow  patterns. 

Figure  11  shows  the  evolution  of  the  z-component  vorticity  u,  in  the  horizontal  plot¬ 
ting  window  (z  =  15  cm)  at  different  time.  After  the  initial  disturbances  have  decayed, 
a  uniform  vorticity  distribution  is  created  at  t  =  120  sec  with  positive  distribution  (solid 
lines)  on  one  side  and  negative  distribution  (dash  lines)  on  the  other.  The  WG  in  the 
corresponding  region,  as  shewn  in  Figure  12,  interacts  with  the  vorticity  field  and  cre¬ 
ates  different  vorticity  tendencies.  The  vorticity  is  enhanced  or  diminished  depending  on 
whether  it  is  in-  or  out-of-phase  with  the  WG  as  shown  in  Figure  13.  The  WG  at 
t  =  120  sec,  being  positive  in  the  vicinity  of  the  center  of  the  wake  column,  interacts  with 
the  positive  vorticity  (solid  lines  in  Figure  11)  on  one  side  of  the  center  line  to  produce  a 
vorticity  stretching  that  enhances  the  positive  vorticity  on  that  side.  The  same  positive 
WG  also  interacts  with  the  negative  vorticity  (dashed  lines  in  Figure  11)  on  the  other  side 
of  the  center  line  to  generate  a  vorticity  stretching  that  reinforces  the  negative  vorticity 
distribution  on  the  other  side.  In  other  words,  the  positive  WG  at  t  =  120  stc  reinforces 
both  the  positive  and  negative  vorticity  strength. 

When  the  WG  becomes  negative,  it  interacts  with  the  wake  vorticity  and  generates  a 
vorticity  stretching  that  enhances  or  diminishes  the  vorticity  strength  depending  whether 
the  stretching  is  out-of-  or  in-phase  with  the  positive  (or  negative)  vorticity. 

The  wake  column  with  the  z-component  vorticity  distribution  as  shown  in  Figure 
11,  though  non-propagating,  is  transported  very  slowly  by  the  induced  current.  As  time 
lapses,  the  initial  vorticity  field  is  further  modulated  by  the  WG,  forming  isolated  vorticity 
centers.  The  modulation  transforms  the  uniform  vorticity  distribution  at  t  =  120  sec 
in  Figure  11  to  two  rows  of  parallel  and  aligned  isolated  (relatively)  vortices  as  shown 
schematically  in  Figure  14a.  As  demonstrated  by  von  Kaman  in  his  origin  study  of  the 
Kaman  vortex  streets,  these  two  rows  of  aligned  and  isolated  vortices  are  unstable  and  tend 
to  gain  stability  by  moving  into  an  alternative  pattern  (Figure  14b)  as  in  the  Kaman  vortex 
streets.  This  re-alignment  leads  to  a  sinusoidal  flow  pattern  (as  observed  in  laboratory 
experiments)  at  t  =  360  stc  as  shown  in  the  vorticity  distribution  in  Figure  11.  in  the 
velocity  field  in  Figure  4,  in  the  trace  concentration  in  Figure  5.  and  in  the  corresponding 
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surface  signatures  in  Figures  7  &  8.  As  the  sinusoidal  patterns  develop,  the  VVG  is  twisted 
away  from  its  original  vertical  direction,  and  re-aligns  with  the  concentrated  vorticity  at 
t  =  360  sec  as  shown  in  Figure  12.  This  re-aligned  VVG  then  continues  to  interact  with  and 
to  strengthen  the  local  vorticity  as  shown  in  the  distribution  of  the  vorticity  stretching  at 
the  same  time  in  Figure  13.  The  present  vorticity  stretching  has  changed  from  its  uniformly 
positive  and  negative  distribution  on  both  sides  of  the  centerline  at  t  =  120  sec,  to  two 
rows  of  aligned  but  relatively  isolated  distributions  at  t  =  360  sec.  This  concentrated 
distribution  has  alternated  positive  and  negative  strength,  and  its  interaction  with  the 
vorticity  is  becoming  different  than  that  at  t  =  120  sec.  As  shown  in  Figure  13  the 
vorticity  stretching  at  t  =  360  sec  has  the  same  positive  or  negative  sign  on  both  sides 
of  the  center  line  while  the  vorticity  distribution  at  t  =  360  sec  is  still  mostly  positive  on 
one  6ide  and  negative  on  the  other,  with  some  sinusoidal  oscillations  superimposed.  This 
vorticity  stretching,  contrast  to  its  earlier  role  at  t  =  120  sec,  strengthens  the  vorticity 
on  one  side  and  weakens  the  vorticity  on  the  other  side  of  the  center  line  at  a  certain 
i  location,  further  encouraging  the  development  of  the  isolated  vortex  pattern  along  the 
x-direction.  This  process  is  pushing  the  flow  to  evolve  towards  a  pattern  similar  to  that  of 
the  Kaxnan  vortex  streets. 

As  the  concentrated  vortex  patterns  further  re-align,  the  VVG  and  the  induced  vortic¬ 
ity  stretching  will  follow  the  re-alignment.  The  induced  vorticity  stretching,  also  alterna¬ 
tively  in  pattern  as  shown  in  Figure  13  for  t  =  600  and  840  sec,  continues  to  amplify  both 
the  positive  and  negative  vorticity  on  both  sides  of  the  centerline  and  further  reinforces 
the  development  of  the  horizontal  vortices.  The  corresponding  flow  patterns  as  shown  in 
the  tracer  and  velocity  fields  then  reveal  the  patterns  reminiscent  of  the  Kaman  vortex 
streets  as  observed  in  experiments. 

In  addition  to  the  organized  patterns  revealed  by  the  vorticity  distribution,  the  mag¬ 
nitude  of  the  vorticity  contours  also  reflects  the  formation  of  the  coherent  structures  in 
the  stratified  late  wake.  The  magnitude  of  the  vorticity  contours  ranges  from  -0  072  to 
0.064/sec  at  120  sec,  -0.04  to  0.035/sec  at  360  sec,  -0  032  to  0.032/sec  at  600,840,  and 
1080  sec,  -0.032  to  0.028/sec  at  1320  and  1560  sec.  and  -0.03C  to  0.028/sec  at  1S00  sec. 
The  change  of  magnitude  between  t  =  120  and  t  =  300  sec  are  quite  distinct  and  reflects 
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the  wake  structure  change  before  and  after  the  onset  of  the  vortex  development.  By  the 
time  the  sinusoidal  patterns  begin  to  form,  the  vorticity  strength  show  little  change  as  an 
indication  of  the  existence  and  persistence  of  the  horizontal  vortices  within  the  flow  field. 

The  proposed  mechanism  for  the  evolution  of  the  horizontal  vortices  is  based  on  the 
vorticity  stretching  in  the  vertical  direction.  This  stretching  is  a  result  of  the  interaction 
of  local  vorticity  produced  by  the  wake  remnants  of  a  submerged  moving  body,  and  the 
vertical  velocity  gradient  resulting  from  the  buoyance  effect  of  the  density  stratification. 
The  interaction  becomes  the  driving  force  that  generates  and  reinforces  the  evolution  of 
the  alternating  and  concentrated  vortex  patterns  as  observed  in  laboratory  experiments. 
This  mechanism  explains  the  delay  appearance  of  the  horizontal  vortices  in  late  wakes, 
and  their  vital  dependency  on  stratified  fluids. 

6.  Summary 

The  organized  flow  structure  generated  by  a  submerged  object  moving  in  a  stratified 
fluid  is  simulated  by  a  three-dimensional  time-dependent  numerical  model.  The  flow  is 
assumed  to  be  incompressible  and  hydrostatic  with  the  Boussinesq  approximation  for  a 
weakly  stratified  fluid.  The  horizontal  turbulent  mixing  is  modeled  by  the  Smagorinsky  for¬ 
mula  to  admit  the  rate  of  deformation,  whereas  the  vertical  mixing  is  based  on  a  Richardson 
number  closure  to  account  for  the  density  and  shear  effects  in  stratified  fluids.  Preliminary 
results  for  the  velocity,  density,  and  tracer  fields  consistently  show  the  existence  and  evolu¬ 
tion  of  such  organized  flow  structures  within  a  stratified  flow  field.  The  surface  signatures 
consequently  reveal  the  development  of  such  vortex  structure  underneath. 

A  mechanism  based  on  the  interaction  of  the  wake  vorticity  and  the  buoyancy  effect 
in  a  stratified  fluid  is  proposed  for  the  generation  and  growth  of  the  horizontal  vortices 
in  stratified  fluids.  The  resultant  vorticity  stretching  is  the  driving  force  for  the  evolution 
of  the  organized  patterns,  from  meandering  to  the  fully  rolled-up  structures  reminiscent 
of  the  Kaman  vortex  streets.  A  test  case  for  homogeneous  fluids  showed  no  development 
of  horizontal  vortices.  This  mechanism  explains  the  delayed  appearance  of  the  horizontal 
vortices  and  their  vital  dependency  on  stratified  fluids.  Furthermore,  our  results  suggest 
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that  simple  two-dimensional  models  cannot  adequately  describe  the  motion  and  character¬ 
istics  of  the  horizontal  vortices.  Only  models  in  which  the  vertical  stratification  is  properly 
represented  are  suitable  for  simulating  these  vortices. 

In  this  paper  only  one  case  of  a  series  of  numerical  experiments  has  been  presented 
and  discussed.  Various  sensitivity  tests  can  be  conducted  to  isolate  effects  of  different 
parameters  such  as  the  depth  of  the  fluid,  the  strength  of  the  stratification,  the  size,  the 
shape,  the  speed  and  the  depth  of  the  moving  object,  etc.  These  and  other  results  will  be 
presented  in  an  upcoming  paper. 
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Figure  4.  Continued. 
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Figure  6.  Density  anomaly  in  a  horizontal  plotting  window  (j-y  plane  at  z 
the  middle  of  the  channel. 
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Figure  6.  Continued. 
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Figure  12.  Continued. 
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Figure  14.  Transition  of  two  rows  of  aligned  and  isolated  vortices  to  gam  stability. 


